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Abstract 

Transmission through a potential barrier, and the related issue of particle 
production from a parametric resonance, are topics of considerable general in- 
terest in quantum physics. The authors have developed a rather general bound 
on quantum transmission probabilities, and recently applied it to bounding the 
greybody factors of a Schwarzschild black hole. In the current article we take a 
different tack — we use the Miller-Good transformation (which maps an initial 
Schrodinger equation to a final Schrodinger equation for a different potential) 
to significantly generalize the previous bound. 
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1 Introduction 



Consider the Schrodinger equation, 



u(x)" + k(x) 2 u(x) = 0, 



(1) 



where k(x) 2 = 2m[E — V(x)]/h 2 . As long as V(x) tends to finite (possibly different) 
constants V±oo on left and right infinity, then for E > max{V +OQ , V-oo} one can set 
up a one- dimensional scattering problem in a completely standard manner — see for 
example [HEIEIIUEIEIITIE]. The scattering problem is completely characterized 
by the transmission and reflection amplitudes (t and r), though the most important 
aspects of the physics can be extracted from the transmission and reflection probabil- 
ities (T = \t\ 2 and R = \r\ 2 ). Relatively little work has gone into providing general 
analytic bounds on the transmission probabilities, (as opposed to approximate es- 
timates), and the only known result as far as we have been able to determine is 
this: 

Theorem 1. Consider the Schrodinger equation (Q]j. Let h(x) > be some positive 
but otherwise arbitrary once-differentiable function. Then the transmission probabil- 
ity is bounded from below by 



To obtain useful information, one should choose asymptotic conditions on the func- 
tion h(x) so that the integral converges — otherwise one obtains the true but trivial 
result T > sech 2 oo = 0. (There is of course a related bound in the reflection prob- 
ability, R, and if one works with the formally equivalent problem of parametric 
oscillations, a bound on the resulting Bolgoliubov coefficients and particle produc- 
tion.) 

This quite remarkable bound was first derived in [5J, with further discussion and 
an alternate proof being provided in [10J. These bounds were originally used as a 
technical step when studying a specific model for sonoluminescence [TTJ, and since 
then have also been used to place limits on particle production in analogue space- 
times [12] and resonant cavities [13], to investigate qubit master equations [TJ], and 
to motivate further general investigations of one-dimensional scattering theory [15J. 
Most recently, these bounds have also been applied to the greybody factors of a 
Schwarzschild black hole [IB] . 



T > sech 2 




y/(/l') 2 + (g - h 2 ) 2 

2h 



dx 



(2) 
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A slightly weaker, but much more tractable, form of the bound can be obtained 
by applying the triangle inequality. For h(x) > 0: 



T > sech 



ln(/i)'| + 



\k 2 -h 2 



dx > . 



(3) 



Five important special cases are: 

• If we take h = k^, where k^ = lim a .^ ±00 k(x), then we have 0, [10 



T > sech^ 



1 



2k 



+oo 



I ^oo ^ I dx 



oo J— CO 



(4) 



If we define k +QO = \im x ^ +ocl k(x) ^ k_ oa = lim, r ^_ 00 k(x), and take h(x) to be 
any function that smoothly and monotonically interpolates between k-oo and 
k +OQ , then we have 



T > sech^ 



In 



boo I t,2 



k 2 -h 2 



dx > . 



(5) 



This is already more general than the most closely related result presented 
in P HO]. 



If we have a single extremum in h(x) then 



T > sech 



In 



h h 

h 2 

"-ext 



-oo |;„2 



h 2 



h 



dx > . 



(6) 



This is already more general than the most closely related result presented 

in pun]. 

If we have a single minimum in k 2 (x), and choose h 2 = max{k 2 , A 2 }, assuming 
£^ in < A 2 < /c 2 ^, (but still permitting fc^ in < 0, so we are allowing for the 
possibility of a classically forbidden region), then 



T > sech < - In 



1 , ( k-i-nr,k- 



A 2 



2A ' U ' |d ''>' 



(7) 



A 2 >k 2 



This is already more general than the most closely related result presented 
in PdOj. 
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• If k 2 (x) has a single minimum and < fc min < fc^, then 

r ,^{i>n(^)}. (S) 

This is the limit of (GO) above as A — > /c min > 0, and is one of the special cases 
considered in [9]. 

In the current article we shall not be seeking to apply the general bound (j2j), its 
weakened form ([3]), or any of its specializations as given in (jlJ)-(jHJ) above. Instead 
we shall be seeking to extend and generalize the bound to make it more powerful. 
The tool we shall use to do this is the Miller-Good transformation [T7] . 



2 The Miller— Good transformation 

Consider the Schrodinger equation and consider the substitution [T7] 

u{x) = -yl== U(X(x)). (9) 
v {%) 

We will want X to be our "new" position variable, so X(x) has to be an invertible 
function, which implies (via, for instance, the inverse function theorem) that we need 
dX/dx 0. In fact, since it is convenient to arrange things so that the variables 
X and x both agree as to which direction is left or right, we can without loss of 
generality assert dX/dx > 0, whence also dx/dX > 0. 
Now compute (using the notation U x = dU / dX) : 

u\x) = U X {X) Vx 1 - \-^j- 2 U{X), (10) 

and 

u"{x) = u xx (x) (x>r - u + \^fr 2 u. (11) 

Insert this into the original Schrodinger equation, u{x)" + k(x) 2 u(x) = 0, to see that 

{k 2 1 X'" 3 (X") 2 1 

irF-2OT + iWr =9, (12) 

which we can write as 

U xx + K 2 U = 0, (13) 
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with 

That is, a Schrodinger equation in terms of u(x) and k(x) has been transformed into 
a completely equivalent Schrodinger equation in terms of U(X) and K(X). You can 
also rewrite this as 

K^^{e + ^(£j'}. (15) 

The combination 

v^f 1 V- 1X "' + 3(X " )2 f!6) 

shows up in numerous a priori unrelated branches of physics and is sometimes re- 
ferred to as the "Schwartzian derivative" . 

• As previously mentioned, to make sure the coordinate transformation x <-> X 
is well defined we want to have X\x) > 0, let us call this j(x) = X\x) with 
j(x) > 0. We can then write 

Let us suppose that ]im x _>± O0 j(x) = j± OQ ^ 0; then K± OQ = k± OCJ /j± OQ , so if 
k 2 (x) has nice asymptotic behaviour allowing one to define a scattering prob- 
lem, then so does K 2 (x). 

• Another possibly more useful substitution (based on what we saw with the 
Schwartzian derivative) is to set J(x)~ 2 = X'(x) with J(x) > 0. We can then 
write 

X 2 = J 4 j^ + yj (18) 

Let us suppose that \im x ^± OQ J(x) = J± OQ ^ 0; then K ±OQ = k± 00 J± 00 , so if 
k 2 (x) has nice asymptotic behaviour allowing one to define a scattering prob- 
lem, so does K 2 (x). 

These observations about the behaviour at spatial infinity lead immediately and 
naturally to the result: 

Theorem 2. Suppose j±oo = 1, (equivalently, J± OQ = 1). Then the "potentials" 
k 2 (x) and K 2 (X) have the same reflection and transmission amplitudes, and same 
reflection and transmission probabilities. 
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This is automatic since K±oo = k± OQ) so equation (QQ) and the transformed equation 
(fl3l) both have the same asymptotic plane-wave solutions. Furthermore the Miller- 
Good transformation ([9]) maps any linear combination of solutions of equation (CQ) 
into the same linear combination of solutions of the transformed equation ([TBI . QED. 

Theorem 3. Suppose j±oo 7^ 1, (equivalently, J±oo 7^ 1). What is the relation 
between the reflection and transmission amplitudes, and reflection and transmission 
probabilities of the two "potentials" k 2 (x) and K 2 (X)? This is also trivial — the 
"potentials" k 2 (x) and K 2 (X) have the same reflection and transmission amplitudes, 
and same reflection and transmission probabilities. 

The only thing that now changes is that the properly normalized asymptotic states 
are distinct 

exptyk^x) exp^K^x) 



but map into each other under the Miller-Good transformation. QED. 



(19) 



3 Improved general bounds 

We already know 

r I-+CO 

T > sech 2 \ ddx\ . (20) 



Here T is the transmission probability, and d is the function 



vW±EZM! r?-n 

with h(x) > 0. But since the scattering problems defined by k(x) and K(X) have 
the same transmission probabilities, we also have 

T>sech 2 |y ^dxj, (22) 

with 

dX = X' dx = j dx, (23) 
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and 



i) 



y/(h X ) 2 + [K 2 - /i 2 ] 2 



2h 



1 



h! 



2h V V X 
1 



r 



2 j 4 y 



h 2 



2hj 



J I 2 j 4 j 2 



(fc') 2 + 



(24) 
(25) 

(26) 



That is: V7i(x) > 0, Vj(x) > we now have (the first form of) the improved bound 

-OO 



T > sech 



J I 2 j 4 j 2 



dx > . (27) 



Since this new bound contains two freely specifiable functions it is definitely stronger 
than the result we started from, (j2J). The result is perhaps a little more manageable 
if we work in terms of J instead of j . We follow the previous logic but now set 



dX = X' dx = J- 2 dx, 



(28) 



and 



(29) 



That is: \/h(x) > 0, VJ(a;) > we have (the second form of) the improved bound 

-oo 



T > sech^ 



1 

2fc 



(h>) 2 + 



J 2 lk 2 + 



J" 



h 2 



J J 2 



dx 



(30) 



A useful further modification is to substitute h = HJ 2 , then WH(x) > 0, VJ(x) > 
we have (the third form of) the improved bound 



T > sech^ 



1 
2H 



H' + 2H 



7 



+ 



j" 

k 2 + H 2 

J 



~1 2 



dx 



(31) 



Equations (l2"7j) . (1501) . and (15T1) . are completely equivalent versions of our new bound. 
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4 Some applications and special cases 



We can now use these improved general bounds, (127)1 . (1301 . and (I3TI) . to obtain 
several more specialized bounds that are applicable in more specific situations. 



4.1 Schwartzian bound 

First, take h = (constant) in equation (130|) . then 



T > sech^ 



A; 2 + 



J 



J2 



dx 



(32) 



In order for this bound to convey nontrivial information we need lim^-i-oo J k = h , 
otherwise the integral diverges and the bound trivializes to T > 0. The further 
specialization of this result reported in P [10] and equation (Tjfl) above corresponds to 
J = (constant) = a/V^oo, which clearly is a weaker bound than that reported here. 
In the present situation we can without loss of generality set h —>■ k^ in which case 



T > sech^ <^ - 
2 



J 



J 2 



dx 



(33) 



We now need lim^-too J = 1 in order to make the integral converge. If k 2 > 0, so 
that there is no classically forbidden region, then we can choose J = yk^Jk, in 
which case 



T > sech 2 <^ - 



1 

7^ 



i 

7^ 



dx 



(34) 



This is a particularly elegant bound in terms of the Schwartzian derivative, [equa- 
tion fj!6p ], which however unfortunately fails if there is a classically forbidden region. 
This bound is also computationally awkward to evaluate for specific potentials. Fur- 
thermore, in the current context there does not seem to be any efficient or especially 
edifying way of choosing J(x) in the forbidden region, and while the bound in equa- 
tion (133]) is explicit it is not particularly useful. 



4.2 Low-energy improvement 

We could alternatively set H = (constant) in equation (I3TI) . to derive 



T > sech 




7 



AH 2 



J" 

k 2 + — -H 2 

J 



dx 



(35) 
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In order for this bound to convey nontrivial information we need lim^-too k 2 = k 2 ^ = 
H 2 , lim^ioo J' = 0, and lim^-too J' = 0. Otherwise the integral diverges and the 
bound trivializes to T > 0. Thus 

T > sech 2 



rj 




2 1 \ 




7_ 


+ [ 



J" 

/-' 2 + -j - kl 



dx 



(36) 



Again, the further specialization of this result reported in [HI QUI and equation (j3J 
above corresponds to J = (constant), which clearly is a weaker bound than that 
reported here. To turn this into something a little more explicit, since J(x) > we 
can without any loss of generality write 



J(x) = exp 



x(x) dx 



where xi x ) is unconstrained. This permits is to write 



T > sech 

Then by the triangle inequality 
T > sech 2 



Ak 2 



[k 2 + X 2 - X' ~ k lof dx 



\X\ 



2,k n 



\k 2 + x 2 



X 



dx 



A further application of the triangle inequality yields 



T > sech 



X + 



2k r 



2kr 



k 2 



x 2 -kU 



dx 



(37) 



(38) 



(39) 



(40) 



Now if k 2 < k 2 ^, (this is not that rare an occurrence, in a non-relativistic quantum 
scattering setting, where k 2 ^ — k 2 = 2mV/h 2 and we have normalized to Voo = 0, it 
corresponds to scattering from a potential that is everywhere positive), then we can 
choose x 2 — k 2 ^ — k 2 so that 



T > sech 



\X\ + 



2k r 



\X 



dx 



(41) 



Assuming a unique maximum for x (again not unreasonable, this corresponds to a 
single hump potential) this implies 



T > sech z 



+ 



lo - k 2 dx 



(42) 
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This is a new and nontrivial bound, which in quantum physics language, where 
k 2 = 2m(E — V)/h 2 , corresponds to 




T > sech 2 {\/^+ / dx } . (43) 

If under the same hypotheses we choose x = 0, then the bound reported in [9l [TO] 
and equation PJ above corresponds to 

T> -^ 2 ^L^^- (44) 



Thus for sufficiently small E the new bound in equation ( 1431) is more stringent than 
the old bound in equation (jSJ) provided 



r— 1 f°° V2mV 



(45) 



Note the long chain of inequalities leading to these results — this suggests that these 
final inequalities fj4*2|) and (|43p are not optimal and that one might still be able to 
strengthen them considerably. 

4.3 WKB-like bound 

Another option is to return to equation (j4*0|) and make the choice y 2 = max{0, —k 2 } = 
k 2 , so that k = \k\ in the classically forbidden region k 2 < 0, while k = in the 
classicallty allowed region k 2 > 0. But then equation PU]) reduces to 

:/'>S(<dr<J / Kdx + ^ + ^+ [ \ k2 °°~ fe2[ d.r } . (Hi) 

fc 2 <0 fc 2 >0 

Key points here are the presence of f k2<0 K dx, the barrier penetration integral that 
normally shows up in the standard WKB approximation to barrier penetration, K max 
the height of the barrier, and L the width of the barrier. These is also a contribution 
from the classically allowed region (as in general there must be, potentials with no 
classically forbidden region still generically have nontrivial scattering). Compare this 
with the standard WKB estimate: 

' wkb ~ sech 2 { J k dx + In 2 } . (47) 

k 2 <0 



10 



This form of the WKB approximation for barrier penetration is derived, for instance, 
in Bohm's classic textbook [18] . and can also be found in many other places. Under 
the usual conditions applying to the WKB approximation for barrier penetration we 
have J k 2 <0 K dx 3> 1, in which case one obtains the more well-known version 

TwKB~expJ -2 J Kdx\. (48) 

I k 2 <0 ) 

The bound in equation (Hoi) is the closest we have so far been able to get to obtaining 
a rigorous bound that somewhat resembles the standard WKB estimate. Again we 
do not expect the bound in equation (146!) to be optimal, and are continuing to search 
for improvements on this WKB-like bound. 



4.4 Further transforming the bound 



In an attempt to strengthen the inequalities (142]) and (143]) . we again use the fact that 
J(x) > to (without any loss of generality) write J(x) = exp [ J x( x ) dx] , where 
x(x) is unconstrained. The general bound in equation ( ]3T1) can then be transformed 
to: For all H(x) > 0, for all xi x ) : 



T > sech z 



OO J 
-rvi 2 



— + 2y 
H 



2 , [k 2 + X 2 + X' ~ gf , 
+ Jp 



(49) 



This leaves us with considerable freedom. Regardless of the sign of k 2 (x), we can 
always choose to enforce k 2 + x 2 — H 2 = 0, and so eliminate either x or H, obtaining 



T > sech' 



OO -j^ 



— + 2VH 2 - k 2 
H 



+ 



{VH 2 - k 2 )' 



H 2 



dx 



(50) 



(subject to H(x) > and H 2 (x) — k 2 (x) > 0), and 



T > sech z 




+ 2x 



, (x') 2 dx 

X 2 + k 2 



(51) 



(subject to x 2 ( x ) + k 2 (x) > 0), respectively. Finding an explicit bound is now largely 
a matter of art rather than method. For example if we take 



H 2 



max{A; 2 ,A 2 } or 



X 



max{0, A 2 - k 2 } 



(52) 
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then from either equation ( J50l) or equation (IBTl) . again under the restriction that we 
are dealing with a single-hump positive potential, we obtain 

T > sech 2 J -In (^p) + (v^!EEW + | v^3^2 dx I . (53) 

Note that A is a free parameter which could in principle be chosen to optimize the 
bound, however the resulting integral equation is too messy to be of any practical 
interest. This bound is somewhat similar to that reported in equations ([7j) and fj42l) . 
but there are some very real differences. 



5 Summary and Discussion 



The bounds presented in this note are generally not "WKB-like" — apart from the 
one case reported in equation ( H6|) there is no need (nor does it seem useful) to 
separate the region of integration into classically allowed and classically forbidden 
regions. In fact it is far from clear how closely these bounds might ultimately be 
related to WKB estimates of the transmission probabilities, and this is an issue to 
which we hope to return in the future. 

We should mention that if one works with the formally equivalent problem of a 
parametric oscillator in the time domain then the relevant differential equation is 

u(t) + k(t) 2 u(t) = 0, (54) 

and instead of asking questions about transmission amplitudes and probabilities one 
is naturally driven to ask formally equivalent questions about Bogoliubov coefficients 
and particle production. The key translation step is to realize that there is an 
equivalence [9l [10] : 

1 1 - T 
T <-> ; N <-> . (55) 

This leads to bounds on the number of particles produced that are of the form 
N > sinh 2 {(some appropriate integral)}. 

To be more explicit about this our new improved bound can be written in any of 
three equivalent forms: 

• For all H(x) > 0, for all J(x) > 0, 



T > sech 




H' + 2H 



7 



-1 2 



+ 



j" 

k 2 + — -H 2 

J 



dx 



(56) 
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For all h(x) > 0, for all J(x) > 0, 



T > sech^ 



2h 



{h'f + 



J 2 I k 2 



J" 
1 



h 2 
J 2 



1 2 



dx 



(57) 



For all h(x) > 0, for all j(x) > 0, 



T > sech^ 



[-00 1 



J l 2 j 4 j 2 ' 



dx } . (58) 



The equivalent statements about particle production are: 
• For all H(t) > 0, for all J(t) > 0, 



N < sinh 2 



OO J 



H' + 2H 



1 



+ 



A; 2 



T 



n 2 



if 2 



df 



(59) 



For all h(t) > 0, for all J(t) > 0, 



N < sinh 2 



2/i 



{h') 2 + 



j 2 U 2 + 



J" 1 h 2 



J 



1 2 



J 2 



dt 



(60) 



For all fc(t) > 0, for all j(t) > 0, 



N < sinh 2 



+oo 



2h 



j I 2 j 4 j 2 < 



(/i') 2 + 



dt } . (61) 



In closing, we reiterate that these general bounds reported in equations (1271) . (1501) . 
and (EH), their specializations in equations ([33]), (J5JJ, (S2J), (S3]), (Pj), and (1351). 
and the equivalent particle production bounds in equations (l5"9l- (j6"Tl) . are all general 
purpose tools that are applicable to a wide variety of physical situations [TT] [12], [13] 
HI [15] [TO] . Furthermore we strongly suspect that further generalizations of these 
bounds are still possible. 
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